Nuclear Multifragmentation Critical Exponents 



In a recent Letter [^ the EoS collaboration presented data of fragmenta- 
tion of 1 A GeV gold nuclei incident on carbon. By analyzing moments of 
Q>^ ! the fragment charge distribution, the authors claim to determine the values 
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G^ I of the critical exponents 7, /?, and r for finite nuclei. These data represent 

a crucial step forward in our understanding of the physics of nuclear frag- 
mentation. However, as we will show in the following, the analysis presented 
^ . in |l[] is not sufficient to support the claim that the critical exponents for 

nuclear fragmentation have been unambiguously determined. 

The main difficulty in observing critical behavior in nuclear fragmentation 
is that nuclei cannot be prepared and held near the temperature and density 
^ ■ associated with the critical point. Instead, complicated nuclear reactions 

^ ■ must be used to excite the nuclei, which may then expand to conditions 

sufficiently close to the critical point. The time spent at these conditions 
during the reaction is an open question. However, even if one assumes that 
conditions sufficiently close to critical are explored in some of the observed 
CK ' reactions, there remain at least two problems with interpreting the subse- 

quent experimental signals. One is to identify which of the resulting parti- 
cles have participated in the equilibrated system near the critical conditions, 
O I and which have resulted from the "pre-equilibrium" stage of the reaction. 

^ I The other problem is to measure the "temperature" of the decaying system. 

These two problems are interrelated in the procedure used in Ref. [|ll, where 



^ 



X 



the authors assume that the observed multiplicity can be used as an indicator 
of temperature. In this comment, we use the percolation model of nuclear 
5^ I fragmentation |0 to demonstrate the nature of these problems in determining 

critical exponents. 

In percolation models one uses a bond breaking parameter p with values 
between and 1, including pc , the "critical value". In this model, near the 
critical point, the charge of the largest fragment is .^max oc (p — Pc)'^, and the 
second moment of the mass distribution is M2 (x \p — Pc\^"', with /3 = 0.41 
and 7 = 1.8. 

We follow the analysis of |l[] and use rric = 26, in accordance with the 
cut employed in [^. It is worth noting that we also find roughly identical 
values of 7 and 7' for the liquid and gas branches, if we examine lnM2 vs. 
ln|m — md and use rric = 26. However, this value of rric = 26 is lower than 
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m{pc), and the numerically extracted value of 7 = 7' is approximately 1.2 
to 1.4, significantly lower than 1.8, the critical exponent for the percolation 
model. These two observations already indicate possible problems in trying 
to find the critical exponents via the analysis employed in |jl|. 

One would only expect ^max oc (m — rric)'^ and M2 oc |m — md'^, if 
multiplicity and p are related in a strictly linear way. Percolation calculations 
show, however, that on average m rises monotonically with p, but not linearly 
near the critical value. Furthermore, for a given value of p, there is always 
a distribution in the values of m. Thus, any translation of an observable, 
O, in terms oi p into terms of m involves a non-trivial convolution, 0{m) = 
J dpm{p) ® 0{p). 

The results in Fig. 1 illustrate the difficulty in extracting the critical 
exponent (3. Here we plot InZ^ax vs. ln|r?T, — md- The open circles represent 
the result from a percolation model with Z = 79 charges. 10^ events were 
generated to provide statistics comparable to the experimental data. The 
best fit to these points results in a slope '/9' = 0.55. For comparison, the 
solid line has a slope of /? = 0.41, the nominal value for percolation. This 
demonstrates that the slope extracted from the present plot is not the critical 
exponent j3. 

In addition to this problem, one expects that a portion of the observed 
multiplicity is comprised of pre-equilibrium particles, which are not part of 
the equilibrated system. In cascade simulations, the number of such charges 
varies from to possibly 15 for the system used in Ref. ^. If it is not pos- 
sible to experimentally (via momentum space analysis, for example) remove 
these, then complicated modelling of the pre-equilibrium reaction dynamics 
is needed. To illustrate the infiuence of such pre-equilibrium contamination 
on the extraction of meaningful values for the critical exponent j3 we show the 
results of a calculation in which we assume that 10 pre-equilibrium particles 
exist along with an equilibrated system of 69 charges. The remaining system 
of 69 charges is then fragmented, again by using the percolation model. Since 
the total observed multiplicity includes the pre-equilibrium particles, plot- 
ting InZmax vs. ln|m — md results in the crosses in Fig. 1. For comparison, 
the dashed fine has a slope of '/?' = 0.29. (Coincidentally, this is the value 
extracted from the experimental data.) 

Similarly, when following the steps of analysis in |]1| for determination 
of 7, we find significant contamination of the 'liquid' branch due to pre- 
equilibrium emission. In particular, we find that the extremely large offset in 



ln|M2| between the liquid and gas branch (about 3.5) in Fig. 2 of |1| seems 
to be caused by this contamination, in combination with the procedure of 
dropping the largest fragment on the "liquid" side of the critical point and 
keeping it on the "gas" side. 

Since the analysis of [|l| does not work for a simple model with known 
values of critical exponents, one should not expect it to yield the correct 
critical exponents for the data, either. 

In summary, while we applaud the effort and beautiful data of |jl| , we do 
not agree with the conclusion that the critical indices of nuclear fragmenta- 
tion have been determined yet. 
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Fig. 1. Percolation model simulation of the fragmentation of gold. Circles: 
Z = 79 system fragmenting; crosses: Z = 69 system fragmenting plus 10 
pre-equilibrium protons. For comparison, the solid line has a slope oi P = 
0.41, and the dashed line '/?' = 0.29. 



